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We study theoretically the exchange of angular momentum between electromagnetic and elec- 
trostatic waves in a plasma, due to the stimulated Raman and Brillouin backscattering processes. 
Angular momentum states for plasmon and phonon fields are introduced for the first time. We 
demonstrate that these states can be excited by nonlinear wave mixing, associated with the scat- 
tering processes. This could be relevant for plasma diagnostics, both in laboratory and in space. 
Nonlinearly coupled paraxial equations and instability growth rates are derived. 
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It is well known that the angular momentum of elec- 
tromagnetic radiation contains two distinct parts, one 
associated with its polarization state, or photon spin, 
the other being the external or orbital photon angular 
momentum (0AM). With the advent of laser beams, an 
increasing interest is being given to the study of photon 
0AM, and various optical experimental configurations 
have been considered [TJ H] El S]. It is now well un- 
derstood that collimated electromagnetic beams, such as 
laser or radio wave beams, can be described by Laguerre- 
Gaussian functions, which provide a natural orthonormal 
basis for a generic beam representation. Utilization of 
photon 0AM states in the low frequency (< GHz) ra- 
dio wave domain was also recently proposed in Ref. [Sj 
as a new method for studying and characterizing radio 
sources in astrophysics. 

The possibility of remote study of space plasma vor- 
ticity by measuring the 0AM of radio beams interacting 
with vortical plasmas was pointed out in Ref. El and a 
more detailed theoretical analysis was given recently by 
studying the electromagnetic wave scattering from the 
plasma medium, with the associated 0AM exchanges be- 
tween the plasma and probing photon beams [7] . A more 
speculative work was also recently published where the 
strong similarities between photon and neutrino disper- 
sion relations were explored, and 0AM states of neutrino 
beams interacting with dense plasmas were considered 

Here we consider the important problem of stimulated 
Raman and Brillouin backscattering of collimated elec- 
tromagnetic beams with finite 0AM in a plasma. This 
also leads us to consider, to our knowledge for the first 



'Electronic address: titomend@ist.utl.pt 

t Also at LOIS Space Centre, Vaxjo University, SE-351 95 Vaxjo, 
Sweden 



time, the possible existence of plasmon and phonon states 
with finite orbital angular momentum. Raman and Bril- 
louin scattering instabilities are well known in the con- 
text of laser fusion [9^ , as possible sources of anomalous 
plasma refiectivity. Raman backscattering is now recog- 
nized as a dominant process for ultra-intense laser plasma 
interactions, in the context of inertial fusion research |10j . 
In all these studies, angular momentum in general, and 
photon 0AM in particular, have been systematically ig- 
nored. On the other hand, there seems to be experimen- 
tal evidence of 0AM dependence in Brillouin scattering 
of radio waves in the ionosphere [TT] , which awaits for a 
deeper theoretical understanding. 

In contrast with the traditional theoretical approach 
[9], we consider the case of an electromagnetic pump 
beam with finite transverse dimension, and arbitrary 
0AM states are considered. Our formalism is quite sim- 
ple but general, and includes the nonlinear coupling be- 
tween incident and backscattered waves with electron 
plasma waves (for the Raman instability) or ion acoustic 
waves (for the Brillouin scattering). Nonlinear parax- 
ial equations for the electromagnetic wave modes will be 
coupled to paraxial equations for the electrostatic wave 
modes. These later equations will allow us to introduce 
the angular momentum states for both plasmons and 
phonons. Notice that the spin effects are absent for the 
electrostatic oscillations, which means that (in contrast 
with the photon case) their 0AM coincides with their 
total angular momentum. 

In the following, instability growth rates for generic 
0AM states of the incident or collimated pump beams 
will be determined, for an infinite and homogeneous 
plasma. An important point to notice is that, even in 
the infinite plasma, our theoretical model will predict a 
well localized backscattering instability, taking place at 
the focal region of the incident beam. 

In isotropic and homogeneous plasmas, where the ions 
are assumed immobile, transverse and longitudinal wave 
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coupling can be described by the following wave equa- 
tions 



^ c^V^ + Ljt]A = — —nA 



pe 



(1) 



and 



where A = A{f, t) is the vector potential describing the 
transverse electromagnetic waves, n the electron density 
perturbations associated with the electrostatic waves, uq 
the equilibrium electron density, UJpe — (e^no/eorTi)^/^ is 
the electron plasma frequency, — e and m are the electron 
charge and mass, cq the vacuum permittivity, c the speed 
of light, and Sg = \/3Te/m the electron thermal velocity, 
for a plasma with electron temperature T^. 

We now assume wave propagation along the z- 
direction, by considering wave solutions of the form 
A = J2j=i 2 ^3 &'^p{ikjZ — iujjt) + C.C., where ujj and kj 
are the frequencies and wavenumbers of the two elec- 
tromagnetic wave modes (the incident and the scattered 
one). Similarly, we can use for the electrostatic oscilla- 
tions h = fiQ exp{ik'z — iuj't) + c.c. We assume that the 
wave amplitudes Aj and no are slowly varying on space 
and time scales much longer than the respective wave- 
lengths and periods. Using such an assumption, we can 
reduce the wave equations ([T|) and ^ to three coupled 
equations of the form 

DiAi = Lul^{no/no)A2 , £12^42 = cj2^(nS/"o)^i (3) 
D'no ^ no(e2fc'V2m2) (^Ai ■ A'^^ 

where we have used the following operators 
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for J = 1, 2, we have assumed the energy and momentum 
conservation relations loi = LO2+L0' and ki ~ k2 + k' , and 
that the linear dispersion relations are satisfied, for the 
transverse and electrostatic modes, fc|c^ = {ojj — OJ^^), 

and k"^Sl = (w'^ — Wpg), respectively. We should keep 
in mind that uj' ~ LUpe, and consequently the incident 
wave frequency is larger than twice this value, loi > 2ujpe. 
It is known that, in order to make all these conditions 
compatible with each other, we have to assume that the 
two electromagnetic wave modes propagate in opposite 
directions, with fci > for the incident wave, and k2 — 
— 1^21 < for the backscattered wave, with fc' > for the 
electrostatic wave. Let us first discuss equations ^ in 
the linear approximation, where the coupling terms on 
the r.h.s. are neglected. The temporal dependence of 
the amplitudes Aj and no will vanish, and these linear 
equations will reduce to pure paraxial equations, taking 




FIG. 1: Amplitude of the electron density fluctuations asso- 
ciated with plasmon states of finite (orbital) angular momen- 
tum, for I' = 1 and p' — 0, in arbitrary units. 



the form 
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Using cylindrical coordinates f = {r,ip,z), the parax- 
ial wave solutions for the first of these equations can 
be written as linear combinations of modes Aj(r) = 
Ap-^i^{z)Fp^^i.{r,z)e''-^'^, where Fp^^i.{r,z) are Laguerre- 
Gaussian functions, with integers Pj,lj representing the 
radial and azimuthal (quantum) numbers, as defined by 
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where Wj = Wj{z) are the beam waists of the 
two electromagnetic wave modes, j = 1,2, Xj = 
{r/wjY, and Lp{x) are associated Laguerre poly- 
nomials. Due to the well known orthogonal- 
ity relations of the Laguerre-Gaussian modes, we 
can write /J" rdrFp.^i, (r, z)F;^ i^_ (r, z) " d^pe'^'^-'-^^ = 
^pi,P2^hM^ where the deltas represent Kronecker sym- 
bols. It is also known that these solutions of the vec- 
tor potential Aj correspond to well defined photon or- 
bital angular momentum states, characterized by the az- 
imuthal quantum numbers Ij . Similarly, we can say that 
the no solutions correspond to a superposition of plas- 
mon angular momentum states, characterized by the ra- 
dial and azimuthal quantum numbers p' ,1' , of the form 
Fpi,ii{r, z)e'' The beam waist w' = w'{z) of the corre- 
sponding plasmon modes is considered here to be of the 
same order as the beam waist of the dominant electro- 
magnetic wave mode, w' Wi. Notice however that the 
electrostatic waves carry no intrinsic angular momentum, 
because in contrast to the transverse photons, plasmons 
have spin zero. This means that the plasmon angular 
momentum states coincide with their total angular mo- 
mentum states. The electric field perturbations associ- 
ated with these plasmon states will remain purely elec- 
trostatic, if they satisfy V x i? = 0. For solutions of the 
form E(r, t) — Fir) exjp(il'ip -\- ik' z — iu't), this will im- 
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ply the existence of angular and radial field components 
such that = I'Ez/kr, and Er = —idEz/k'dr, where 
Ez is the axial field component. The same will be true 
for the electron mean velocity v — —i(eQUj'/enQ)E. In the 
present calculations we only need to use the explicit form 
for the electron density fiuctuations ho, as described by 
the paraxial equation ([5|. The corresponding solution 
can be described as a superposition of orthogonal LG 
modes 



il' ifi ik' z — iu' t 



Cz + C.C. (7) 



These plasmon modes can then take quite unusual 
shapes, as illustrated in figure 1. Going back to the cou- 
pled nonlinear equations (3), we can see that they de- 
scribe an exchange of energy, linear momentum, as well 
as orbital angular momentum, between the two electro- 
magnetic waves and the longitudinal waves. Replacing 
the linear solutions (with amplitudes now depending on 
time) in these equations, and integrating over the radial 
variable r, we arrive at the following nonlinear coupled 
equations 



dai / dt = —iCin' d2 , dSydt = iC2n'dl 
dn'/dt = ~iC'{di ■ dl) 

with the nonlinear coupling coeflicients 



(8) 
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We have also used the simplified notation dj — Ap. ^i. (2, t) 
and n' = npi^ii{z,t), and introduced the quantity R = 
R{z) = /p Fp^^i^Fp^^i^Fpi ^iirdr. In order to understand 
the physical meaning of these coupled equations for the 
Laguerre-Gaussian modes, let us consider the important 
case of stimulated Raman scattering, by considering an 
intense incident wave with amplitude Oi , and by using the 
parametric approximation dai/dt ~ 0. We also assume 
maximum coupling conditions, corresponding to parallel 
polarization di \\ a2. From the two remaining equations, 
we can then easily derive 
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This leads to the following unstable solutions n'{z,t) — 
n'(z, 0) exp(7t), and 02(2, i) = 02(2, 0) cxp(7t), with 
growth rate 7. They show that the incident wave with 
amplitude ai can excite electron plasma waves, or plas- 
mons, with different angular momentum states. In the 
case of plasmons carrying no angular momentum, the 
backscattered wave will be in a state I2 = —h- The ob- 
servation of a backscattered signal with I2 ^ —h will 
therefore reveal the existence of plasmons with non-zero 
angular momentum states propagating in the medium. 
Raman backscattering can then be used as a powerful di- 
agnostic method to detect internal plasma vorticity. Let 
us now consider the interesting situation where well de- 
fined plasmon angular momentum states can be excited 
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FIG. 2: Normalized amplitude of the plasmon mode excited 
across the laser focal region, for three different times. 



from the outside. This can be done by using two counter- 
propagating electromagnetic waves. Going back to equa- 
tions we can determine the amplitude of the excited 
plasmon state, as a function of the initial backscattered 
wave signal 02(2,0). Again, we assume |ai| ^ |a2| and 
use the parametric approximation. Equations ([s]), for 
parallel polarization, can then be reduced to 



da*2 
dt 



2n a^e 



dn' 
~dt 



aia2e (11) 



Here we have introduced a finite frequency mismatch 
Auj = LOi ~ u)2 — Lo' , because the two electromagnetic 
wave modes imposed from the outside are not necessar- 
ily in the optimum matching conditions. We can then 
derive 
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again with 7^ = C"C2|aipi?^. For initial conditions such 
that n'{z, 0) = and 02(2, 0) is arbitrary (but still obey- 
ing the linear paraxial solution along 2), we can write the 
solution of these equations in the form 

02(2, t) = a2{z, 0) cosh(5t)e*'^'^*/2 (^3) 
with <7 = [72 - (Aw/2)2]i/2^ and 

n'(z, t) = -^ — 02(2, 0) sinh((7i)e*^'^*/2 (14) 
9 

These solutions describe the growth of the backscattered 
signal with orbital angular momentum state I2, and the 
excitation of a plasmon angular momentum state char- 
acterized by the azimuthal number /' — h — l2- Notice 
that the above growth rates depend on the axial position 
2. This means that, very soon, the axial profile of both 
the excited backscattered and electrostatic modes will 
deviate from its linear solution, as illustrated in figure 2. 
However, the radial beam profile will remain unchanged. 

Stimulated Brillouin scattering with orbital angular 
momentum can be treated in a similar ground. In this 



4 



case, the nonHnear plasmon equation ([2| is replaced by 
the nonlinear equation 
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where h now represents the electron density oscillations 
associated with the ion acoustic waves, Vs = ^ZT^jM 
is the ion acoustic velocity, and Ze and M are the ion 
charge and mass. Following the same procedure, and 
assuming that n = exjp(ik' z — iuj't) + c.c., we can arrive 
at the nonlinear paraxial equation for the ion acoustic 
wave amplitude 
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where the differential operator is the same as before, but 
with the electron thermal velocity replaced by the 
ion acoustic velocity Vg, and with the parallel wavenum- 
ber k' satisfying now the ion acoustic dispersion relation 
k'vs ~ Lo' . In the linear approximation, we recover the 
electrostatic paraxial equation, similar to the second of 
equations ^ . This shows that angular momentum states 
for ion acoustic waves, formally identical to those of plas- 
mons, can be defined such that nQ{r) — n'Fp' ii{r, z)e*' 
with n' ~ np'^ii{z,t). Inserting these linear solutions, in 
the nonlinear paraxial equation (151, we get for the time 



dependent ion acoustic wave amplitude 
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This leads to stimulated Brillouin backscattering so- 
lutions, with a growth rate 7^, determined by 7^ — 
CBC2|aip. All the qualitative features discussed above 
for the Raman backscattering case can be repeated here. 
In particular, the use of electromagnetic wave scattering 
as a diagnostic probe of internal plasma vorticity, and 
the excitation of well defined angular momentum states 
of the phonon spectrum by two counter-propagating elec- 
tromagnetic waves. 

In this work we have considered stimulated Raman 
and Brillouin backscattering of electromagnetic waves in 



a plasma with photon orbital angular momentum. We 
have also introduced, to our knowledge for the first time, 
the related concepts of plasmon and phonon angular mo- 
mentum states. The field modes associated with these 
electrostatic wave modes are determined by solutions of 
appropriate paraxial equations, similar to those describ- 
ing coUimated electromagnetic wave beams near the focal 
region. In the case of these electrostatic quasi-particles, 
the orbital angular momentum coincides with the to- 
tal angular momentum, because, in contrast to photons, 
plasmons and phonons have no spin. 

We have discussed the nonlinear coupling between in- 
cident and backscattered radiation in uniform plasmas, 
and derived the corresponding growth rates, in the para- 
metric approximation. These results generalize the well 
known results of stimulated Raman and Brillouin scat- 
tering, by including the finite size and radial profile of 
the wave modes, corresponding to the various orbital an- 
gular momentum states. We have shown that an addi- 
tional selection rule for nonlinear wave interaction, asso- 
ciated with the conservation of angular momentum, was 
added to the usual energy and linear momentum selec- 
tion rules. In particular, we have shown that, by using 
two counter-propagating electromagnetic waves with well 
defined orbital angular momentum, we can excite specific 
states of non-zero plasmon and phonon angular momen- 
tum. Experimental verification of our theoretical model 
could lead to a novel result on basic plasma physics, to 
be considered in the future. 

In this work we have only considered three wave cou- 
pling processes. But the present formalism can easily be 
extended to four wave coupling, which could be particu- 
larly interesting for the case of Brillouin scattering. No- 
tice that, although the present model only appUes to uni- 
form plasmas, the unstable region where stimulated scat- 
tering takes place is highly localized, first because of the 
finite transverse width of the electromagnetic beams, sec- 
ond because of the enhancement of the instability growth 
rate in the axial direction, over a distance eventually 
much shorter than the Rayleigh length. 
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